Self-driven nonlinear dynamics in magneto-optical traps 
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We present a theoretical model describing recently observed collective effects in large magneto- 
optically trapped atomic ensembles. Based on a kinetic description we develop an efficient test 
particle method, which in addition to the single atom light pressure accounts for other relevant 
effects such as laser attenuation and forces due to multiply scattered light with position dependent 
absorption cross sections. Our calculations confirm the existence of a dynamical instability and 
provide deeper insights into the observed system dynamics. 
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Since its first realization in 1987 [ij the magneto- 
optical trap (MOT) has become a standard technique 
for providing a robust source of large numbers of cold 
atoms. While multiple scattering of the absorbed laser- 
light is known as a major limitation for achieving Bose- 
Eintein condensation, it also leads to interesting collec- 
tive effects which have been studied over the last years 
0, 0, 0, 0, 0, and a variety of static structures has 
been observed and investigated by different theoretical 
approaches |1,I1)IEI3- 

Only recently, experiments have revealed a so far unex- 
plored dynamical instability in three-dimensional MOTs 
connected with the appearance of self-excited radial os- 
cillations 1 1 1| , which constitutes a complex nonlinear dy- 
namics phenomenon. Understanding the observed effect 
turns out to be of broader interest, as it provides a clean 
laboratory realization of similar plasma- and astrophysi- 
cal phenomena, such as, e.g., pulsating stars 12], which 
are generally difficult to access. 

Here we develop a theoretical model, describing the ob- 
served instability and providing a physical picture of the 
underlying mechanism. As discussed in J\, sub-Doppler 
cooling mechanisms only affect a very small fraction of 
large trapped atom clouds. Hence, the overall behavior 
of large atomic ensembles is well described within a ba- 
sic Doppler-cooling picture, where the resulting trapping 
force along each laser beam can be written as J^] 

Hlpi^^v) = -J- is+(T+ix,v) - s_CT_(r,i;)) , (1) 

where 

..=ao(l + 3(., + ._)+4i^i^^^)" (2) 

is the absorption cross section for the two laser beams 
(including a saturation by the 3 pairs of laser beams), 
CTo = 3A/27r the on-resonance absorption cross section, 
A the laser wavelength, F the transition linewidth, S 



the detuning from resonance, fj,x determines the Zee- 
man shift of the atomic transition due to the MOT 
magnetic field and s± = /±//sat denotes the satura- 
tion parameter of the respective laser beam of inten- 
sity I± with /sat being the saturation intensity of the 
atomic transition. For the discussion below it is con- 
venient to split the force according to Ftrap(2^,i') — 
a{\x\)x + P{\x\,\v\)v, with a{\x\) = i^trap(|a;|, 0)/|a::| and 
/3(|x|, |t;|) = (Ft.ap(N|, \v\) - Ft,MM/\v\- 

In order to simplify our theoretical considerations we 
use the following spherical symmetric generalization of 
eq.((TJ) 

Ftrap = a(r)r + /3(r,'y)v , (3) 

While experimental confinement configurations generally 
do not obey this symmetry, eq.||2Jl describes the impor- 
tant features of the resulting force in both the linear and 
nonlinear trapping regions. 

At higher densities, attenuation of the laser light inside 
the cloud, results in an additional effective confining force 
experienced by the atoms To account for this effect 
within our spherical symmetry assumption, the spatial 
intensity profile is obtained from 

s_ = soe^ '^o" '^+(''')''('~')-/J''^-('^')p('-') (4) 

where sq is the saturation parameter of the incident 
beam. Moreover, multiple scattering of the absorbed 
laser light inside the cloud leads to an additional outward 
directed pressure, caused by an effective interaction be- 
tween the atoms j^. Neglecting higher order scattering 
events, which are known to screen the atom-atom inter- 
action [l^ . a photon scattered off an atom at position 
ri exerts an average force on an absorbing atom at r2 
according to 0, Q 

„ 3/sat / (+) , (-)\ r2 — ri 
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The reabsorption cross section aisc is obtained by con- 
volving the absorption cross section of the emitted Hght 
with the emission spectrum of the atom at ri in the pres- 
ence of either left or right circularly polarized laser light. 
Note that Crsc niay depend on both coordinates via the 
space dependence of the local laser intensities as well as 
of the respective detunings. Previously |3,0,0,S0E|j 
such coordinate dependencies have been neglected, which 



J 



according to eq.© results in a Coulomb-like interaction 
with effective charges, again underlining the close anal- 
ogy with plasma and gravitational physics problems. In 
large clouds, however, we find the position dependence of 
the effective charges to be important for the static and 
dynamics properties of the trapped atom cloud. 

Starting from eqs.(|21)-© the collective system dynam- 
ics is described by the following kinetic equation 
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for the atomic phase space density /(r,v,t), where 

F,„f(r)= y"F,,e(r',r)/(r',v)dvdr (7) 

and M is the mass of the atoms. Heating by sponta- 
neous emission and photon exchange [l6l ll^ has been 
neglected, since for the densities considered in this work 
the corresponding thermal pressure is much smaller than 
the pressure resulting from the effective atomic repulsion. 

Note that eq.Q goes beyond a local-density approxi- 
mation [ifij . retaining the complete position dependence 
of (Trsc and the density dependence of Frsc- In fact, this 
nonlocal space dependence of all forces in eq. © in addi- 
tion to their local dependence on the atom position ren- 
ders a direct numerical solution of eq.® very demanding. 
Alternatively, we apply an efficient numerical procedure 
based on a test-particle treatment, similar to particle- 
in-cell methods [l8j, frequently used for plasma physics 
problems. More specifically, we represent the atomic den- 
sity by an ensemble of Nt < 10^ test particles, whose 
number is typically chosen to be less than the actual par- 
ticle number to reduce the numerical effort. The respec- 
tive absorption cross sections and masses of the test par- 
ticles are adjusted, such that the results are independent 
of the number Nt of test particles. By propagating every 
particle according to the forces eq. Q and O we obtain 
the time dependent density from which we calculate the 
local intensities and the resulting forces to advance the 
next timestep. 

To study its stationary properties we evolve the atomic 
cloud until it relaxes to the selfconsistent, stationary solu- 
tion of eq.©, which we found to exist only below a criti- 
cal atom number iVc. Fig^Ji shows the calculated station- 
ary density profile for N = 1.15 x 10^ Rubidium atoms 
and typical MOT-parameters of / = l.OmW/cm^, S = 
— l.Sr and T/ fj, = 4.7mm (corresponding to 9G/cm)^3|- 
As can be seen, the calculated density is well described 
by a truncated Gaussian profile. As the atom number 
is decreased the truncation radius R decreases relative 
to the rms-width of the corresponding Gaussian, ulti- 



mately leading to a transition into a uniform density pro- 
file. Similar changes in the density profile have also been 
reported in MOTs, where the nonlinearity of the poten- 
tial arises from sub-Doppler trapping mechanisms [tIIi^. 
In the present case the observed transition results from 
the nonlinearity and the position dependence of the re- 
absorption cross section and, hence, can not be found 
under the assumption of linear trap ping forces and pure 
Coulomb-like interactions |3.l3.ll0|. 

Let us now turn to the most striking result of our calcu- 
lations. As we further increase the number N of atoms 
the cloud becomes unstable at a critical atom number 
A^c, corresponding to a critical radius Re- By varying 
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FIG. 1: (a) Radial density profile (solid line) of trapped 
Rubidium atoms with S — — l.SF, / = 1.0m W/cm and 
T/fi = 4.7mm together witfi a Gaussian fit (dashed line), (b) 
Radial dependence of Fmi (dashed line), Ftrap (dotted line) 
and i^mf + i^trap (solid line), (c) Radial dependence of the 
damping constant /3{r,v = 0). The gray shaded area marks 
the region of active atomic motion. 
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FIG. 2: Relaxation of the MOT's RMS-radius after switching 
the detuning from Sq = — 1.55r to 5 = — l.SF for two different 
particle numbers of iV = 9 x 10* (a) and iV = 1.3 x 10^ 
(b). The remaining parameters are the same as in fig0 The 
solid line in (a) shows a fit of Act sin(cjt + (j})e~^^'^ + CToo to 
the numerical data (circles). The solid line in (b) shows the 
numerical data in the oscillating regime of N > Nc- 
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FIG. 3: (a) Calculated MOT size as a function of the atom 
number (circles) fitted by a power- law dependence (solid line) . 
The dotted line corresponds to the critical radius i?c = 5/fi 
and the dashed line in both figures indicates the critical atom 
number A'c beyond which the dynamical instability sets in. 
(b) Real part of the MOT's stability coefficient (circles) fitted 
by TCx{N - iVc)" with Nc = 1.226 x 10^ and k = 0.55. 



the various MOT-parameters, we find that the critical 
radius is uniquely determined by the relation Rc = 6/ u 
(see figE^), confirming the conclusion reached in 
This fact is illustrated in figCJ) and^, where we show 
the radial dependence of the trapping and interaction 
force as well as the damping constant /3(r, u = 0). The 
damping constant /3(r, 0) reverses its sign at i?c = S/fi. 
Hence, any small velocity of atoms outside of Rc will 
be enhanced. While inward moving particles will be 
damped again when entering the negative-/? region, out- 
ward moving atoms around Rc are further accelerated 
away from the trap center, since the single atom light 
pressure force is largely balanced by the interaction force 
around r = Rc- Their motion around the fixed point 
(r — Rc^v = 0) will become unstable and limited by 
the non linear terms of the force. At larger distances, 
however, the total force reverses its sign again, since the 
interaction force decreases much more rapidly than the 
trapping, due to the radially increasing Zeeman shift (see 
fig. ^p). Hence, if during the expansion, the atoms did 
not acquire a velocity beyond the capture range of the 
MOT, a stable limit cycle will be reached. 

In order to characterize the onset of the instability 
we analyze the clouds RMS radius a = \J (r^) /3 and 
study its sensitivity against a small perturbation. More 
precisely, we start from a stationary density correspond- 
ing to some detuning which is instantly increased to 
(5 (closer to resonance), leading to damped oscillations 
of a towards its new equilibrium value <7oo, as shown 
in fig|2^. From a fit to a damped harmonic oscillation 
a = Acre"*/"^ sin(a;t -I- 0) -I- a^o we obtain the damp- 
ing time T and frequency to corresponding to the real 



and imaginary part of the respective Lyapunov exponent 
A = -f iu (figEt)- For increasing N and the param- 
eters of figI21the instability sets in at an atom number of 
Nc = 1.226 X 10^ and with a critical exponent of 0.55. 
On the other hand, the frequency of the cloud oscilla- 
tion evolves continuously through the instability thresh- 
old (see figEJl, indicating that the onset of the instability 
proceeds via an supercritical Hopf-bifurcation. 

A reduction of the system properties to a single quan- 
tity like the cloud's RMS-radius is clearly helpful for un- 
derstanding the transition into the oscillating regime. On 
the other hand the fully resolved space-time evolution of 
the atomic density such as shown in figQl reveals much 
more detailed information about the complicated dynam- 
ics of the cloud. Indeed, the complex density patterns at 
larger atom numbers (see fig0J) shows that the oscillation 
dynamics is much more complex than a simple breathing 
mode, as suggested by the simple size oscillations close 
to the instability threshold (see figOb). 

In fact, the oscillation is triggered by an outer fraction 
of atoms, which gain energy as they move in and out of 
the active region of r > _Rc, which is indicated by the 
horizontal white line in fig^ When bouncing back on 
the low-energetic atoms, the gained energy is deposited 
by exciting a density wave just inside the region with 
/3(r, 0) < 0. Subsequently, the formed nonlinear excita- 
tion propagates towards the trap center along the diag- 
onal blue line drawn in fig^and thereby loosing energy, 
mostly due to the damping by the cooling lasers. As can 
be seen in fig 01 this not only leads to a flattening and 
broadening of the density wave until it disappears, but 
also to a deceleration as indicated by the deviation of 
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optical confined atoms provide an ideal laboratory sys- 
tem for further exploration of the rich spectrum of self- 
driven motion, including variable system geometries, ef- 
fects of external driving and possibilities to control the 
system dynamics. 
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FIG. 4: (color online) Spatio-temporal evolution of the atomic 
density for A'^ = 1.7 x 10®. The MOT parameters are the same 
as in figQ The lines are discussed in the text. 



the moving maximum from the blue line at smaller dis- 
tances. At the same time, the edge region of the atomic 
cloud starts to relax, causing some atoms to be again 
accelerated away from the center and the whole process 
repeats itself. Although this scenario clearly provides 
the basic mechanism for the observed oscillations, our 
calculations reveal a number of finer details (see figEJ 
still to be understood. Moreover, additional damping 
mechanisms, similar to Landau-damping of plasma waves 
might also play a role for the system dynamics, raising 
the interesting question of how the present nonlocal po- 
sition dependence of the effective charges manifests itself 
in known plasma kinetic effects. 

In conclusion, large clouds of magneto-optical confined 
atoms have been found to exhibit a very complex nonlin- 
ear dynamics. Our theoretical description has revealed 
the onset of a deterministic instability connected with 
self-sustained oscillations in agreement with recent exper- 
iments jOl . It has been found that a number of different 
effects, such as the attenuation of the trap lasers, rescat- 
tering of the absorbed laser light as well as the position 
dependence of the respective absorption cross sections 
are all necessary to explain the observed phenomenon. 
A stability analysis of the MOT size has shown that the 
transition proceeds via a supercritical Hopf-bifurcation. 
The obtained density evolution revealed the build-up of 
complex nonlinear excitations driven by the combined ac- 
tion of the light-pressure force and the effective atomic 
interaction, which results in an active atomic motion at 
large distances. Similar types of active or self-driven mo- 
tion are currently discussed in a broad range of differ- 
ent applications, such as collective swarm dvnamics |2ri. 
propagation of waves or dissipative solitons j22l| in 
reaction-diffusion systems or grain motion in dusty plas- 
mas l^^l- Hence, we believe that large clouds of magneto- 
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